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Abstract: The rapid economic development that the Hotan Oasis in Xinjiang Uygur Autonomous Region, 
China has undergone in recent years may face some challenges in its ecological environment. Therefore, 
an analysis of the spatiotemporal changes in ecological environment of the Hotan Oasis is important for 
its sustainable development. First, we constructed an improved remote sensing-based ecological index 
(RSET) in 1990, 1995, 2000, 2005, 2010, 2015 and 2020 on the Google Earth Engine (GEE) platform and 
implemented change detection for their spatial distribution. Second, we performed a spatial 
autocorrelation analysis on RSEI distribution map and used land-use and land-cover change (LUCC) data 
to analyze the reasons of RSEI changes. Finally, we investigated the applicability of improved RSEI to 
arid area. The results showed that mean of RSEI rose from 0.41 to 0.50, showing a slight upward trend. 
During the 30-a period, 2.66% of the regions improved significantly, 10.74% improved moderately and 
32.21% improved slightly, respectively. The global Moran's I were 0.891, 0.889, 0.847 and 0.777 for 1990, 
2000, 2010 and 2020, respectively, and the local indicators of spatial autocorrelation (LISA) distribution 
map showed that the high-high cluster was mainly distributed in the central part of the Hotan Oasis, and 
the low-low cluster was mainly distributed in the outer edge of the oasis. RSEI at the periphery of the 
oasis changes from low to high with time, with the fragmentation of RSEI distribution within the oasis 
increasing. Its distribution and changes are predominantly driven by anthropologic factors, including the 
expansion of artificial oasis into the desert, the replacement of desert ecosystems by farmland ecosystems, 
and the increase in the distribution of impervious surfaces. The improved RSEI can reflect the 
eco-environmental quality effectively of the oasis in arid area with relatively high applicability. The high 
efficiency exhibited with this approach makes it convenient for rapid, high frequency and macroscopic 
monitoring of eco-environmental quality in study area. 
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1 Introduction 


An oasis is an important part of arid area, and its ecology is relatively independent and fragile 
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(Liu et al., 2018). As a relatively independent geographical unit, an oasis is the carrier of 
anthropological socioeconomic activities in arid areas and has a positive regulatory effect on local 
climate, water resources and ecosystem (Hao et al., 2019). Hotan Oasis is a piedmont oasis with 
limited natural resources such as water, fragile ecology, prominent environmental constraints and 
resource consumption problems in a typical arid area, Xinjiang Uygur Autonomous Region, 
China (Wang et al., 2020). Additionally, the proportion of agricultural land has increased 
substantially since 2010 (Dong et al., 2019), causing the edge of the oasis to advance considerably 
toward the desert. The extent of expansion can be estimated visually from satellite images. 
Balancing socioeconomic development and ecological environment is key to the sustainability of 
oasis development. Therefore, monitoring the distribution and changes of ecological environment 
of the oasis in a timely and rapid manner is particularly important. 

A combination of remote sensing spatial information technology and geoinformation 
technology for spatial observation and analysis, which has the unique advantages of update 
quickly and wide range, is widely used in ecological environment monitoring (Yu et al., 2014). 
Furthermore, many researchers have developed spectral algorithms or indices that are closely 
related to ecological environment. For example, the normalized difference vegetation index 
(NDVI) is widely used in various ecological studies (Cao et al., 2018; Chu et al., 2019; Hassan et 
al., 2019; Fung and Siu, 2000). Leaf area index (LAT) is used to reflect the number of plant leaves, 
canopy structure changes, vitality of plant communities and environmental effects (Fang et al., 
2019; Kganyago et al., 2020). Land surface temperature (LST) is used to study the effects of the 
ecological environment and temperature, such as the impact of urban heat islands on an urban 
ecological environment (Alexander, 2020; Zhao et al., 2020). Other indicators used to 
characterize ecological environment include enhanced vegetation index (EVI) to correct soil 
background and atmospheric noise aspects (Wang et al., 2003), hyperspectral flower index (HFD 
to assess flower cover in the flowering state (Chen et al., 2009), normalized difference water 
index (NDWI; Magruder et al., 2019) and modified normalized difference water index (MNDWI; 
Xu, 2005) that reflect the moisture status of ecological environment. The abovementioned indices 
are all refined single indicators for a certain field with a complex ecological environment. The use 
of single indicators to represent the overall change in the quality of ecological environment would 
be insufficient. 

Many multifactor indicators based on remote sensing have been proposed to represent complex 
ecological environment. As ecological environment is monitored using multiple indicators from 
different aspects, the use of multiple factors is more advantageous than the application of single 
factor (Xu et al., 2019). For example, ecological environmental status index (EI), which uses 
biological abundance, vegetation coverage, water network density and other indicators to 
comprehensively evaluate and systematically analyze ecological environment, is widely used to 
evaluate ecological environment-carrying capacity (Zhu et al., 2017). The synthesized drought 
index (SDI; Li et al., 2020), scaled drought condition index (SDCI; Rhee et al., 2010) and 
improved optimal scaled drought condition index (OSDCI; Guo et al., 2019) are used for 
large-scale drought monitoring using precipitation, temperature and vegetation health conditions 
as input variables. The moderate resolution imaging spectrometer (MODIS) global disturbance 
index (MGDI), which uses LST and EVI as input variables, is also used to monitor the 
disturbance of forest ecosystems (Mildrexler et al., 2009). However, although these multifactor 
indices can reflect the quality of ecological environment in the study area, obtaining each of these 
variables may be challenging, and the weight of each factor needs to be determined repeatedly by 
experiments. In addition, when applied to different study areas, the weights differ (Xu et al., 
2019). In the remote sensing-based ecological index (RSEI) proposed by Xu (2013), he calculated 
four single indicators (greenness, wetness, dryness and temperature, which are closely related to 
ecological environment) based on land satellite imagery and coupled into a comprehensive index 
by principal component analysis (PCA). This model has been proven, many times, to be credible 
and reliable at reflecting the quality of regional ecological environment. For example, results 
showed that the overall wetland eco-environmental quality is closely related to the proportion of 
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humidity in RSEI (Jing et al., 2020). RSEI was used to predict the impact of potential increase in 
population and impervious surfaces on ecological environment (Xu et al. 2018), and discuss the 
spatiotemporal evolution of ecological quality of island under different urban densities (Liu et al., 
2021). 

The aforementioned research on RSEI usually focuses on wet monsoon areas or densely 
populated cities; however, the Hotan Oasis is a typical oasis in an arid area. Ecological 
environment is characterized by high potential evaporation and relatively few water resources. 
Therefore, we attempted to make the following improvements to explore the applicability of RSEI 
to arid area. As farmland is the mainly land type in the Hotan Oasis, we kept the following two 
indicators: greenness and humidity. The proportion of urban impervious surface distribution in the 
Hotan Oasis is relatively small; however, it includes a large expanse of bare land and desert. For 
the dryness index, we removed the original index-based building index (IBI) from RSEI and 
retained the bare soil index (BSI), which is sensitive to bare land and desert. The salinity index is 
more sensitive to reflect the degree of salinization of dry land surface, therefore we replaced 
original heat index with it. Thus, the aforementioned remote sensing indices of greenness, 
humidity, dryness and salinity, which reflect the macro eco-environmental quality of arid area, are 
re-combined via PCA to form a new RSEI for arid area. 


RSEI change 
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Spatial 
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RSEI analysis 
based on LUCC 


Applicability 
verification 


RSEI maps in 1990,1995, 
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Fig. 1 Workflow of the study. GEE, Google Earth Engine; NDVI, normalized difference vegetation index; WET, 
wetness; BSI, bare soil index; SI, salinity index; PCA, principal component analysis; RSEI, remote sensing-based 
ecological index; LUCC, land-use and land-cover change. 
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However, the use of local computer resources to process the aforementioned data would 
become increasingly time consuming as the study area increased. This generates many 
intermediate data processing products that occupy local storage resources. The efficiency of data 
analysis is greatly limited and does not allow for timely follow-up monitoring and analysis of 
ecological environment. 

Google Earth Engine (GEE) is a cloud computing platform that can batch process remote 
sensing images and other big data from the earth observations. Owing to its powerful computing 
capability, GEE can rapidly process multitemporal and large-scale remote sensing data compared 
with traditional desktop remote sensing image processing software (Mutanga and Kumar, 2019). 

We combined the improved RSEI and the GEE platform for the analysis of the temporal and 
spatial changes in the eco-environmental quality of the Hotan Oasis. The workflow is shown in 
Figure 1: (1) we used the dataset of Landsat TM, ETM and OLI/TIRS on the GEE platform to 
construct RSEI; (2) we analyzed the temporal and spatial changes in the eco-environmental quality 
of the Hotan Oasis in 1990, 2000, 2010 and 2020; (3) analysis of the spatial correlation of RSEI in 
different years and revealed the changing trend of its eco-environmental quality; and (4) the 
distribution of RSEI and the reasons for variation were analyzed using LUCC data. The JavaScript 
code used to build RSEI in this study is available in the GEE platform, as follows: 
https://code.earthengine.google.com/?accept_repo=users/changhorn4616/RSEI. 


2 Materials and methods 


2.1 Study area 


Hotan Oasis is located at 36°30’N—38°15'N and 79°05'E-80°40’E, and the terrain is dominated by 
sloping plains with an altitude ranging from 1300 to 1450 m a.s.l. The Yurungkax River and 
Karakax River pass through the oasis and meet at the lower reaches of the two rivers (Guo et al., 
2016). It includes Hotan City, Hotan County, Moyu County, Lop County and Kunyu City. The 
study area has an arid desert climate, and is windy and sandy with little precipitation and high 
potential evaporation (Ding et al., 2020). 

Vast deserts and Gobi surround the Hotan Oasis. These land surfaces have almost no vegetation 
cover, and have not changed significantly over the years. The monitoring and analysis of these 
deserts are not significant for this study, but can have a huge impact on the overall PCA results. 
Therefore, a general boundary was drawn between the Hotan Oasis and the surrounding desert 
through visual interpretation and removed the desert area, leaving a study area of approximately 
5000 km?, as shown in Figure 2. 


2.2 Data source and data processing 


In 2021, U.S. Geological Survey (USGS) launched a new Landsat Collection 2 Level 2 dataset on 
the GEE platform (https://developers.google.com/earth-engine/datasets/catalog/landsat). This 
dataset contains Landsat surface reflectance and surface temperature bands after orthorectification 
(Cook et al., 2014). To better reflect ecological environment conditions of vegetation growing 
season in the study area, we selected images with less than 20% cloud cover from July to August in 
the above dataset for seven years: 1990, 1995, 2000, 2005, 2010, 2015 and 2020. Using the cloud 
removal algorithm provided on the GEE platform, we used the pixel quality band QA_PIXEL to 
remove clouds and their shadows, and used the median value to synthesize the multilayer image 
into a good-quality multiband image. To ensure enough quantity and quality of the images, we 
combined three years of images into one year; for example, the image of 1990 is a combination of 
all images from July to August in 1989, 1990 and 1991 with less than 20% of cloud cover. Owing 
to the large time scale, the satellites/sensors used in different years and datasets varied. As shown 
in Table 1, Landsat 5 TM images were selected for 1990, 1995 and 2010, Landsat 7 ETM+ 
images were selected for 2000 and 2005, and Landsat 8 OLI/TIRS images were selected for 2015 
and 2020. 
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Table 1 Basic information of Landsat data sources 


Nene Parasi Satellite Used band peany seen iste pe Preeti 
1990 LANDSAT/LT05/C02/ Landsat5 SR_B (1, 2, 3,4, 5, 4 1989/1990/1991 146/033 Jul-Aug 
1995 TI_L2 /TM 7) ST_B6 14 1994/1995/1996 146/034 <20% 
2000 LANDSAT/LE07/C02/ Landsat7 SR_B (1, 2,3,4,5, 12 1999/2000/2001 146/035 30m 

2005 TI_L2 ETM+ 7) ST_B6 19 2004/2005/2006 147/034 

2010 "e ~ > ae 345) 13 2009/2010/2011 147/035 

2015 LANDSAT/LC08/C02/ Landsat 8 SR_B(2,3,4,5,6, 12  2014/2015/2016 

2020 TI_L2 JOLITIRS 7) ST_B1O 13. 2019/2020/2021 


We generated LUCC data using unsupervised classification based on the characteristics of the 
Hotan Oasis. Socio-economic data were obtained from the Xinjiang Statistical Yearbook during 
1993-2020 and the China County Statistical Yearbook during 1993-2020. 


2.3 Methodology 


2.3.1 Construction of RSEI on the GEE platform 

This study explored the applicability of the improved RSEI in the Hotan Oasis based on the GEE 
platform to evaluate the eco-environmental quality of the study area. This index is a composite of 
four indicators related to the ecological environment of arid area: greenness, wetness, dryness and 
salinity, which obtained by remote sensing data, and then coupled via PCA; the index is expressed 
as follows: 


RSEI=f (greenness, wetness, dryness, salinity). (1) 
(1) Greenness 
npviz 5 —* , (2) 
N+R 


where N is the near-infrared band of the Landsat data; and R is the red band of the Landsat data. 
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(2) Wetness 

The wetness index reflects the water content of soil, vegetation and artificial ground. The 
tasseled cap transformation algorithm can extract corresponding components from an image (Joshi 
et al., 2019). In this study, the moisture component of tassel cap transformation as a wetness 
indicator to reflect the humidity of land surface in the study area. Images from different sensors had 
different coefficients for each band when calculating tassel cap transformation humidity component, 
which is expressed as follows (Huang et al., 2002; Baig et al., 2014): 


WET(TM)=0.03 1 5ppiue +0202 IPgreen 0-3 102p req 


+0.1594px 7p — 0.6806pswr1 — 0.6109 swiss (3) 
WET(ETM+)=0.2626pp jue +0.2141P reen +0.0926P 
+0.0656pyip — 0.7629 swig; — 0-5388PswR2> (4) 


WET (OLI)=0.1511 Pp jye+0-1973P green +0-3283P pe 


+0.3407 pip — 9.711 7pgwip, — 9-4559Pswiro> () 


where WET(TM), WET(ETM+) and WET(OLI) are the moisture component of tassel cap 
transformation of the Landsat TM, ETM+ and OLI images, respectively, and Ppiue, P green, Pred, PNIR, 
pswiri and pswir2 represent the surface reflectance of blue, green, red, near infrared, shortwave 
infrared 1 and shortwave infrared 2 bands of the Landsat TM, ETM+ and OLI images, 
respectively. 

The presence of surface water affected the load size of ground moisture in the PCA. We selected 
seasonal water and permanent water data from the JRC Yearly Water Classification History dataset 
(https://developers.google.com/earth-engine/datasets/catalog/JRC_GSW1_3_YearlyHistory?hl=en) 
(Pekel et al., 2016) on the GEE platform as mask files and used these to mask the water bodies of 
original images. 

(3) Dryness 

BSI reflects the dryness of ground surface, and is an important indicator that reflects the 
ecological condition of arid area. Therefore, dryness is expressed using BSI, which is expressed as 
follows: 


[(Pswri + Pied ) — (PNR + Prine) 


BSI= . 
[(Pswri t+Prea )+(nnr + Prue) 


(6) 


(4) Salinity 

Salinity is an important factor that limits ecological environment of the oasis in arid area. The 
retrieval of soil salinity indicators from images can reflect soil salinity status effectively (Whitney 
et al., 2018; Yu et al., 2010). Salinity index (SI) inversion equation used in this paper is expressed as 


follows (Khan et al., 2005): 
SEy Potue * Pred: (7) 
(5) PCA 


Four indicators including NDVI, WET, BSI and SI were combined into a comprehensive 
indicator through principal component transformation, and most of the information was 
concentrated on the first two principal components: PC1 and PC2 (Zheng et al., 2020). This way, 
eco-environmental quality monitoring can be realized by means of remote sensing without 
manually configuring the weights of various indicators. Because the range of values for the four 
indicators is not unified, the weight of each indicator becomes unbalanced if PCA is conducted 
directly on the indicators. Therefore, before PCA, the range of each indicator value must be unified 
to 0-1. The normalization formula is expressed as follows: 

I EE I min 
NI, a (8) 


max min 


where NZ, is the index normalization result; J; is the pixel value of each pixel; /maxis the maximum 
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pixel value of the entire image; and Jmin is the minimum pixel value of the entire image. 

The four single-band images obtained after the above normalization were synthesized into a 
multiband image, and the same operation was performed for every year. Seven multiband images in 
1990, 1995, 2000, 2005, 2010, 2015 and 2020 were obtained. If the contribution rates of the 
eigenvalues of PC1 obtained after performing PCA exceed 80% and the eigenvectors are stable, 
PC1 can be considered to integrate most of features of these images. Generally, the larger the value 
of PC1, the more inferior the eco-environmental quality that may be reflected; this is contrary to 
typical logical judgments. Therefore, the positive and negative values of PC1 may be transposed to 
obtain the initial RSET, expressed by RSEIo as follows: 


RSEI,=1— PCI. (9) 


To facilitate quantitative analysis, we normalized RSEIo to obtain the improved RSEI, which can 
reflect the eco-environmental quality of arid area at a macroscopic level. The folowing formula was 
used: 


RSEI=(RSEI, — RSE] 
where RSEI and RSEI 


2.3.2 Spatial autocorrelation analysis 


YRSElo mag — RSET nin )> (10) 


Omin Omax 


omax represent the minimum and maximum values of RSEIp, respectively. 


Omin 
Spatial autocorrelation can reflect the spatial distribution pattern (Soka and Oden, 1978; Legendre, 
1993) and interdependence among elements in the research field (Getis, 2007). The application of 
spatial autocorrelation analysis in ecological environment, especially in large-scale remote sensing 
ecological environment analysis, is relatively successful (Jing et al., 2020; Xiong et al., 2021). In 
this study, we adopted two commonly used spatial statistical methods, global spatial autocorrelation 
analysis and LISA. Global spatial autocorrelation was used to analyze the overall spatial distribution 
characteristics of RSEI in the research area, which was represented by global Moran's 7 (Moran, 
1948; Yang et al., 2021) via the following formulas: 


fet jes Mi ji HIG F) 
global Moran's J= Z Dida : : ; (11) 


So EE ~ xy’ 
So = 2 2 Wi, (12) 


where x; and x, are the spatial attribute values for positions ¿ and j, respectively; x is the average 


value of the spatial attributes of all elements in the study area; n is the total number of space 
elements; W;,; is the spatial weight; and So is the collection of all spatial weights. The range of global 
Moran's J values is (-1, 1); A global Moran's I>0 means that RSEI has a clustered distribution in 
space. The closer it is to 1, the higher the degree of clustering. A global Moran's /<O means that 
RSEI has a discrete distribution in space. The closer it is to —1, the higher the degree of this discrete 
distribution. A global Moran's /=0 means that the attribute is randomly distributed in the overall 
space (Li et al., 2020). 

LISA can reflect the local correlation of the spatial distribution of the attributes of a research 
object. Even if the global Moran's Z is 0, there may be local cluster phenomenon; therefore, a 
combination of LISA and global Moran's J analysis is necessary. The formula is expressed as 
follows (Anselin, 1995): 


LISA=~ (13) 


Ss j=l, jži 
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n =\2 
Ss? 2a. a (14) 


í n-l1 


where x; and x, are the attribute values of spatial elements i and j, respectively; Wj, is the spatial 


weight between spatial elements i and j; x is the mean value of the corresponding attributes; and 
nis the total number of spatial elements. 

The results of local spatial autocorrelation analysis generated a LISA clustering map where five 
spatial distribution types were presented: high-high (H-H) cluster, low-low (L-L) cluster, 
high-low (H-L) outlier, low-high (L-H) outlier and not significant. H-H indicates the presence of 
high values is surrounded by high values and L-L indicates the presence of low values is 
surrounded by low values. H-L indicates a high value abnormality and L-H indicates a low value 
abnormality. A not significant value indicates that the attribute value is close to a random 
distribution. 

2.3.3 Impact of LUCC on RSEI 

Over the past 30 a, the land cover type of the Hotan Oasis has changed significantly, which directly 
leads to changes in RSEI. To investigate how the distribution and changes in RSEI are reflected in 
LUCC, this study uses Landsat imagery from two time points, 1990 and 2020 (the same data 
sources used to calculate RSEI in 1990 and 2020), to classify their land cover types. 


3 Results 


3.1 PCA of RSEI 


The results of the four indices, i.e., NDVI, WET, BSI and SI, obtained using PCA for each year are 
shown in Table 2. The eigenvalue contribution rates of the PC1 in 1990, 1995, 2000, 2005, 2010, 
2015 and 2020 were 89.63%, 88.67%, 89.59%, 86.94%, 91.52%, 90.60% and 91.59%, respectively. 
In other words, they all exceeded 85.00%, indicating that PC1 integrates most of the features of the 
four indices and can be used as the initial RSEI, that is, RSEIo. Take the PCA reports for 1990 and 
2020 as an example, as shown in Table 3, The rates of contribution from the eigenvalues of PC2, 
PC3 and PC4 were all found to be below 10% and the positive and negative directions are unstable, 
reflecting more noisy features. Therefore, they were disregarded in the analysis. In PC1, NDVI and 
WET eigenvectors for each year are negative, and BSI and SI eigenvectors are positive, which 
appears to contradict the actual ecological environment condition. However, PC1 is not the final 
result, and a positive and negative conversion operation for PC1 using Equation 9 is required. In this 
way, the sign of the eigenvector is consistent with the actual ecological environment condition. 


Table 2 First principal component (PC1) of principal component analysis for RSEI in 1990, 1995, 2000, 2005, 
2010, 2015 and 2020 


PCI 1990 1995 2000 2005 2010 2015 2020 


NDVI -0.48641 -0.48552 0.49201 -0.48959 0.49658 -0.49817 -0.50478 
WET -0.47930 -0.47372 0.47392 -0.46008 +~—- -0.48033. -0.47307 -0.47814 
BSI 0.51779 0.52127 0.51721 0.52479 0.51152 0.51710 0.51375 
SI 0.51534 0.51782 0.51557 0.52271 0.51092 0.51052 0.50263 
Eigenvalue 0.09923 0.09809 0.09906 0.09599 0.10061 0.09903 0.10018 
oe of eigenvalue 264 88.67 89.59 86.94 91.52 90.60 91.59 


Note: RSEI, remote sensing-based ecological index; NDVI, normalized difference vegetation index; WET, wetness; BSI, bare soil index; 
SI, salinity index. 


JOURNAL OF ARID LAND 


Table 3 Principal component (PC) analysis for 1990 and 2020 


Year Indicator PC1 PC2 PC3 PC4 
NDVI —0.48641 0.66914 0.10060 0.55275 
WET —0.47930 —0.70834 —0.20970 0.47387 
BSI 0.51779 0.17157 0.63423 0.54792 

i SI 0.51534 0.14517 —0.73733 0.41194 
Eigenvalue 0.09923 0.00942 0.00186 0.00020 
Percentage of eigenvalue (%) 89.63 8.51 1.68 0.18 
NDVI —0.50478 0.55881 —0.10010 0.65032 
WET —0.47814 —0.80865 —0.17005 0.29755 
BSI 0.51375 —0.18263 0.54236 0.63918 

i 0.50263 —0.02139 —0.81665 0.28283 
Eigenvalue 0.10018 0.00625 0.00283 0.00012 
Percentage of eigenvalues (%) 91.59 5.71 2.59 0.11 


Note: RSEI, remote sensing-based ecological index; NDVI, normalized difference vegetation index; WET, wetness; BSI, bare soil index; 
SI, salinity index. 


3.2 Spatial distribution of RSEI 


We processed RSEI maps in 1990, 1995, 2000, 2005, 2010, 2015 and 2020 to calculate the 
multi-year mean of RSEI of the Hotan Oasis, as shown in Figure 3. The high values of multi-year 
mean RSEI correspond to the middle of the oasis and both sides of the river, except for the urban 
built-up areas. The highest values occurring in the south and southeast of built-up areas in Moyu 
County, in a strip-like distribution. The outer side of the oasis and bare land of the river channel are 
generally distributed with low values and surround the entire oasis. The medium and low values 
were mainly distributed in Kunyu City, the northwestern and northeastern parts of the oasis, and the 
built-up areas of Moyu County, Hotan City and Lop County. 

The value of RSEI was divided into five levels: poor, fair, moderate, good and excellent 
corresponding to 0.0-0.2, 0.2—0.4, 0.4-0.6, 0.6—0.8 and 0.8-1.0, respectively, as shown in Figure 
4. In Kunyu City, the distribution was almost always poor in 1990, 1995 and 2000, and fair, 
moderate and good in 2005 and increasingly every subsequent year; their areas increased to the 
maximum in 2020. Combined with Figure 5, the average values of RSEI in 1990, 1995, 2000, 
2005, 2010, 2015 and 2020 were 0.41, 0.42, 0.42, 0.46, 0.46, 0.50 and 0.50, respectively, showing 
an increasing trend. The poor distribution proportion decreased year by year from 26% in 1990 to 
4% in 2020, and the percentages of fair and moderate distributions of RSEI increased from 23% 
to 27% and from 23% to 36%, respectively. The percentage of good distribution of RSEI 
fluctuates, but is relatively stable. The percentage of excellent distribution increased slightly but 
not significantly. 


3.3 Change detection of RSEI 


To explore the spatial and temporal changes in RSEI, we performed a change detection of RSEI in 
the Hotan Oasis during 1990-2000, 2000-2010, 2010-2020 and 1990-2020, as shown in Table 4. If 
RSEI level does not change over the year, it is indicated by 0, for example, from moderate to 
moderate; if the region grows by one level, it is indicated by 1, for example, from moderate to good; 
and if the level drops, for example, from excellent to good, it is indicated by —1. Similarly, if the 
level rises by 2 or 3 and drops by 2 or 3, it is indicated by 2 or 3 and —2 or —3, respectively. We 
omitted 4 and —4 from the analysis because the area of the region ascending by four levels and 
descending by four levels was too small. We named the levels that changed by 3, 2, 1, 0, —1, —2 and 
-3 as follows: obviously improved (OI), generally improved (GI), slightly improved (SD, invariable 
(IN), slightly deteriorated (SD), generally deteriorated (GD) and obviously deteriorated (OD), 
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Fig.3 Spatial distribution of multi-year mean RSEI (remote sensing-based ecological index) in the Hotan Oasis 


respectively. Results showed that IN accounted for 74.35%, 61.55%, 55.38% and 41.11% during 
1990-2000, 2000-2010, 2010-2020 and 1990-2020, respectively, indicating that approximately 
half of the regional eco-environmental quality has changed. The percentages of OI and OD during 
1990-2000, 2000-2010 and 2010-2020 did not exceed 1.00%, and even throughout the 30-a period 
from 1990 to 2020, the percentage of IO was only 2.66%, while the percentage of OD was only 
0.10%, indicating that the area with changes obviously in eco-environmental quality was very small 
(Table 4). During 1990-2000, 2000-2010, 2010-2020 and 1990-2020, the percentages of GD area 
were 0.28%, 0.22%, 0.93% and 1.36%, whereas those of GI were 0.85%, 3.45%, 4.52% and 10.7%, 
respectively. The increase of GI was significantly higher. Similarly, the percentages of SD were 
11.74%, 8.28%, 13.44% and 11.80%, and those of IS were 12.75%, 26.34%, 25.33% and 32.21%, 
respectively. The differences between SI and SD for the aforementioned periods were 1.01%, 
18.06%, 11.89% and 20.41%, respectively. That is, the amount of spatial and temporal changes in 
RSEI of the Hotan Oasis remained stable from 1990 to 2000, but increased significantly between 
2000-2010 and 2010-2020. 
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Fig. 4 Spatial distribution of RSEI (remote sensing-based ecological index) level in the Hotan Oasis in 1990, 
1995, 2000, 2005, 2010, 2015 and 2020 
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Fig. 5 Percentage of different levels of RSEI (remote sensing-based ecological index) in the Hotan Oasis in 1990, 


1995, 2000, 2005, 2010, 2015 and 2020 
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Table 4 Area and percentage of RSEI level change during 1990-2000, 2000-2010, 2010-2020 and 1990-2020 


Year Index 3 (OD 2 (GD 1(SD 0 (IN) —1 (SD) —2 (GD) —3 (OD) 
1990-2000 Area (km?) 0.88 36.59 547.22 3191.04 504.04 12.05 0.34 
Percentage (%) 0.02 0.85 12.75 74.35 11.74 0.28 0.01 
2000-2010 Area (km?) 6.79 147.91 1130.46 2641.94 355.40 9.47 0.18 
Percentage (%) 0.16 3.45 26.34 61.55 8.28 0.22 0.00 
2010-2020 Area (km?) 13.80 194.02 1087.26 2376.80 576.91 39.75 3:53 
Percentage (%) 0.32 4.52 25.33 55.38 13.44 0.93 0.08 
1990-2020 Area (km?) 114.35 460.84 1382.44 1764.33 506.54 58.32 4.23 
Percentage (%) 2.66 10.74 32.21 41.11 11.80 1.36 0.10 


Note: RSEI, remote sensing-based ecological index. OI, GI, SI, IN, SD, GD and OD represent that the RSEI level of the area is 
obviously improved, generally improved, slightly improved, invariable, slightly deteriorated, generally deteriorated and obviously 
deteriorated, respectively. The number of 3—1, 0 and —1— —3 represent the improved, unchanged and deteriorated levels of RSEI, 
respectively. 


RSEI level changes in the four periods are shown in Figure 6. During 1990-2000, most regions 
did not change significantly, and even the changed regions did not exhibit significant distribution 
characteristics. During 2000-2010, SI and GI distributions of larger areas occurred in the central 
part of Kunyu City, the northern part of Moyu County, the eastern and northeastern parts of Lop 
County, and on both sides of the downstream channels of the two rivers. A ring of GD hotspots 
formed at the periphery of the built-up area of Hotan City, and a considerable area of SD 
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Fig. 6 RSEI (remote sensing-based ecological index) level change during 1990-2000 (a), 2000-2010 (b), 
2010-2020 (c) and 1990-2020 (d). OI, GI, SI, IN, SD, GD and OD represent that RSEI level of the area is 
obviously improved, generally improved, slightly improved, invariable, slightly deteriorated, generally 
deteriorated and obviously deteriorated, respectively. The number of 3—1, 0 and —1— —3 mean the improved, 
unchanged and deteriorated levels of RSEI, respectively. 
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distribution occurred in the eastern and central parts of Lop County. SI and GI were distributed in 
the reclamation area in the northern part of Kunyu City, in the outermost part of the two rivers, and 
in the southeastern part of Lop County during 2010-2020. The distribution of GD and OD in the 
eastern and southern parts of the urban built-up area of Moyu County indicates a significant 
degradation of eco-environmental quality. The distribution of degraded areas in the periphery of the 
urban built-up area of Hotan City has expanded sharply to the west. From RSEI level change in the 
30-a period of 1990-2020 shown in Figure 6, eco-environmental quality near the edge of the oasis 
increases, while eco-environmental quality in the middle of the oasis decreases, and the "three cores 
and one line" macroscopic ecological distribution pattern among Hotan City, Moyu County, Lop 
County and their connecting roads has been formed. 


3.4 Spatial autocorrelation of RSEI 


To analyze the spatial correlation of RSEI, we sampled RSEI distribution in 1990, 2000, 2010 and 
2020 with a 1 kmx1 km grid to generate 5956 sampling points in the Hotan Oasis. The scatter plots 
of global Moran's J of RSEI in 1990, 2000, 2010 and 2020 are shown in Figure 7. Most of scattered 
points were distributed in the first and third quadrants, indicating that eco-environmental quality of 
the study area was positively correlated with its spatial distribution. Global Moran's J in 1990, 2000, 
2010 and 2020 were 0.891, 0.889, 0.847 and 0.777, respectively, showing a high degree of 
aggregated distribution characteristic. 
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Fig. 7 Global Moran's / scatterplot of RSEI (remote sensing-based ecological index) level change in the Hotan 
Oasis in 1990 (a), 2000 (b), 2010 (c) and 2020 (d). H-H, L-L, H-L and L-H are high-high cluster, low-low cluster, 
high-low outlier and low-high outlier, respectively. 
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3.5 Local indicators of spatial autocorrelation of RSEI 


To explore the spatial characteristics of RSEI, we generated LISA clustering maps of RSEI in 1990, 
2000, 2010 and 2020, as shown in Figure 8. Compared with 1990, the extent of NS (not significant) 
expanded in 2000, and it continued to expand in 2010, and reached the most in 2020, especially in 
the built-up area of Moyu County, Hotan City and Lop County. The H-H cluster area is divided into 
four parts: the central reclamation area of Kunyu City, the northwestern reclamation area of Moyu 
County, and the western and eastern reclamation areas of the Yurungkax River. The fragmentation 
of the H-H cluster area is significantly higher in 2020 than in previous years; therefore, the global 
Moran's J is the lowest in 2020. The distribution of the H-L outlier and L-H outlier is very small and 
sporadic, indicating that there are few extreme points of eco-environmental quality in the entire 
oasis. Overall, the areas of H-H cluster are mainly distributed in the middle of the Hotan Oasis, and 
the areas of L-L cluster are mainly distributed in the outer edges of the oasis, reflecting the trend 
that the overall eco-environmental quality decreases from the inner part of the oasis to the edges. 
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Fig. 8 LISA (local indicators of spatial autocorrelation) clustering map of RSEI in the Hotan Oasis in 1990 (a), 
2000 (b), 2010 (c) and 2020 (d) 


3.6 Impact of LUCC on RSEI 


We used the Landsat images of the Hotan Oasis in 1990 and 2020 (the same images as those used to 
calculate RSET) to generate land use maps, as shown in Figure 9. According to its characteristics, we 
classified the Hotan Oasis into six types of land use: high coverage cultivated land, general 
coverage cultivated land, forest, grass and shrubland, impervious surfaces, water bodies and bare 
land. High coverage cultivated land is the land on which high covered crops are grown such as rice 
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and wheat (Fig. 9a). General coverage cultivated land is the land on which general covered crops 
are grown, consisting mainly of fruit trees (such as jujube and walnuts) or fruit trees combined with 
food crops (Fig. 9b). Impervious surfaces are mainly urban built-up areas (Fig. 9d). We divided the 
forest, grassland and shrubland into one category (Fig. 9c). The water bodies include river, ditch, 
reservoir and pond (Fig. 9e). The bare land include desert (Fig. 9f). 


Land-use type 
ms High coverage cultivated land 
m= General coverage cultivated land 
m= Forest, grass and shrubland 
m= Impervious surfaces 
m Water bodies 
= Bare land 


0 20 km 


Fig. 9 Spatial distribution map of LUCC (land-use and land-cover change) in the Hotan Oasis in 2000 and 2020 
(a, b, c, d, e and f are the Google Earth images of six types of typical land surfaces. Boxes are the location of 
these land surfaces) 


Compared with the map of RSEI spatial distribution (Fig. 4), we found that the distribution of the 
bare land and the worst eco-environmental quality areas (Denoted by Poor) almost overlap. 
Compared with 1990, the distribution of this land type in 2020 is already very small, and was 
replaced by cultivated land, which correspondingly showed by a decrease in poor distribution of 
RSEI and an increase in moderate distribution of RSEI. 

In 1990, there was almost no cultivated land at Kunyu City, while in 2020, cultivated land 
occupied large area, indicating that agricultural activities were frequent over the 30-a development 
and that the eco-environmental quality has changed significantly. These results confirmed by the 
increase in reclamation areas on both sides of the Yurungkax and Karakax rivers and in the eastern 
part of Lop County. On the other hand, the area of impervious surface undergone a considerable 
increase in 2020 compared with 1990, indicating that the continuous urbanization of the Hotan 
Oasis. Correspondingly, the fragmentation of the distribution of RSEI increased inside the oasis. 


4 Discussion 


4.1 Validation of indicators for RSEI 


Based on the 5956 sampling points, we projected NDVI, WET and RSEI onto a three-dimensional 
space as one group, and BSI, SI and RSEI onto a three-dimensional space as another group, as 
shown in Figure 10. RSEI increases when NDVI and WET increase, but decreases when BSI and 
SI increase, showing different and strong linear relationships. Therefore, the greenness 
represented by NDVI and the wetness represented by WET were positively correlated with the 
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eco-environmental quality represented by RSEI, which is consistent with the actual situation (Fig. 
10a). The similar relationships of NDVI and WET with RSEI used for the study of ecological 
quality in Er Hai Lake Basin by Xiong et al. (2021). The dryness represented by BSI and salinity 
represented by SI are negatively correlated with the eco-environmental quality represented by 
RSET, in line with the actual situation (Fig. 10b). This is similar to the relationship between the 
surface albedo (albedo), SI and soil salinity ecological index (SSEI) from remote sensing of 
ecological environment in salinized irrigation areas in the Manas River Basin in Xinjiang Uygur 
Autonomous Region (Zhang et al., 2017). 


(b) 


> 2 
ate Gd 


e Correlation between NDVI and WET e Correlation between BSI and SI 
e Correlation between RSEI and NDVI ® Correlation between RSEI and BSI 
e Correlation between RSEI and WET ® Correlation between RSEI and SI 


Fig. 10 Three-dimensional scatter plots of NDVI, WET and RSEI (a); and of BSI, SI and RSEI (b). NDVI, 
normalized difference vegetation index; WET, wetness; RSEI, remote sensing-based ecological index; BSI, bare 
soil index; SI, salinity index; 


The question of whether RSEI we constructed for arid area is more representative than single 
indicators can be analyzed using the Pearson correlation coefficient (Feidas et al., 2014). We 
selected four single indicators and RSEI composite indicator for the correlation coefficient matrix 
calculation from 1990 and 2020, respectively, as shown in Figure 11. The red color represent a 
positive correlation. The darker the red, the stronger the positive correlation. The darker the blue, 
the stronger the negative correlation. Mean is the average correlation, using the year 1990 as an 
example, Meannpv1_1990 = (|0.74|+|-0.95|+|-0.90])/3=0.86. The mean of NDVI, WET, BSI, SI and 
RSEI in 1990 was 0.86, 0.85, 0.93, 0.93 and 0.96, and that in 2020 was 0.90, 0.85, 0.93, 0.91 and 
0.96, respectively. Thus, the mean of NDVI, WET, BSI, SI and RSEI in these two years was 0.88, 


Correlation Correlation 
(a) 1990 coefficient (b) 2020 coefficient 
1.00 1.00 
NDVI | NDVI 0.74 -0.95 -0.90 0.94 | ogo 9 NDVI|NDVI 0.80 -0.98 -0.93 0.96 0.80 
0.60 0.60 
WET WET -0.88 -0.93 0.93 0.40 WET WET -0.88 -0.88 0.92 0.40 
0.20 0.20 
BSI BSI 0.95 —0.99 0.00 BSI BSI 0.92 —0.98 0.00 
—0.20 —0.20 
SI SI —0.99 —0.40 SI SI —0.97 —0.40 
—0.60 —0.60 
RSEI RSEI -0.80  RSEI RSEI -0.80 
—1.00 —1.00 

NDVI WET BSI SI RSEI NDVI WET BSI SI RSEI 


Fig. 11 Correlation matrix of each indicator and RSEI in 1990 (a) and 2020 (b). NDVI, normalized difference 
vegetation index; WET, wetness; BSI, bare soil index; SI, salinity index; RSEI, remote sensing-based ecological 
index. 
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0.85, 0.93, 0.92 and 0.96, respectively. For a single indicator, BSI has the highest average 
correlation between indicators, which reaches 0.93 in 1990 and 2020. For the composite indicator, 
the average correlation coefficient between RSEI and these four indicators for each year exceeds 
0.93. In addition to integrating the main information for each indicator, the improved RSEI is more 
representative than any single indicators, better illustrating the eco-environmental quality of the 
Hotan Oasis. 

Yu et al. (2021) studied the salinized land degradation index (SDI) of eco-environmental quality 
in the Aral Sea Basin. Although SDI was positively correlated with SI and albedo and negatively 
correlated with NDVI and land surface soil moisture index (LSM), positive and negative values of 
PC1 after PCA was not converted in the study. According to the final results, the correlation 
coefficients between combined indicators and single indicators were both greater than those 
between single indicators themselves, which is consistent with the result of current study. 

The climate of the Hotan Oasis is dry with strong potential evaporation (Ding et al., 2020). The 
PC1 loadings of BSI and SI, representing surface dryness and salinity, ranged from 0.50 to 0.52 in 
the PCA (Table 2). Those of NDVI and WET that represented surface greenness and humidity, 
ranged from 0.46 to 0.50, with a slightly lower weighting compared with those of BSI and SI. The 
result is similar to that of the Aral Sea Basin (Wang et al., 2021). The PC1 loadings for the 
greenness and moisture index in the constructed arid remote sensing ecological index (ARSED) are 
also relatively low due to the low vegetation coverage in the area (Jiang et al., 2019). Wang et al. 
(2019) conducted a study of ecological quality assessment of the Manas Lake in Xinjiang Uygur 
Autonomous Region. They found that the PC1 loadings of NDVI and WET were only 0.12 and 0.33, 
respectively. The lower values of these loadings can be attributed to the low coverage of vegetation 
in the area, scare precipitation and high potential evaporation (Wang, 2021). Our result found that 
salinization in the oasis-desert transition zone of the Hotan Oasis led to a slightly higher share of SI 
in the PC1 loadings. Similar results were found on the eco-environmental quality of the 
Gurbantunggut Desert in Xinjiang Uygur Autonomous Region (Jiang et al., 2019) and Dulan 
County in Qinghai Province (Jia et al., 2021). 

Different indicators have different degrees of spectral response to ground objects (Xue et al., 
2017). To highlight their sensitivity, we considered the year 2020 as an example and extracted the 
image pixel corresponding to the value range of 0.5-1.0, and then selected Sentinel 2 images with a 
higher spatial resolution as the effect comparison, as shown in Figure 12. In Figure 12, al and a2 are 
mainly distributed in the periphery of urban built-up areas with a high vegetation coverage, which is 
consistent with that reflected in the standard false color image of Sentinel 2 (Fig. 12a6). The areas 
reflected by b1 and b2 (Fig. 12) are cultivated land on both sides of the Karakax River, which has 
the highest crop coverage owing to the stable water supply, while bare ground on both sides of river, 
desert and artificial surfaces are all clearly shown on Fig. 12b3. In Figure 12, cl reflects the striped 
reclamation area in the northwestern part of the Hotan oasis, bare land in the middle is relatively 
low and flat, and the degree of salinization is higher than elsewhere, as reflected in the c4 of Figure 
12. 


4.2 Validation of spatiotemporal distribution for RSEI 


We selected representative indicators from the Xinjiang Statistical Yearbook and the China County 
Statistical Yearbook from 1992, 1995, 2000, 2005, 2010, 2015 and 2019 to generate the statistical 
table of socio-economic data from the Hotan Oasis, as shown in Table 5. Results showed that 
population, total power of agricultural machinery, cultivated land area and grain production have 
increasing trends. Population growth leads to an increase in the demand for cultivated land (Deng 
et al., 2015). The reclamation of new cultivated land occurs at the edge of the Hotan Oasis, 
transforming desert into cultivated land. This result explains the increasing RSEI values at the 
edge of the Hotan Oasis (Chen et al., 2019). 

The statistical table shows that gross regional product (GRP) for both the secondary industry 
GRP and the tertiary industry GRP are growing rapidly. The growth in secondary and tertiary 
industries, combined with the emergence of large numbers of field roads and settlements in rural 
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Fig. 12 Comparison among indicators. a, b and c represent different sampling locations; 6 and 7 represents the 
standard false color image and true color image of Sentinel 2, respectively. NDVI, normalized difference 
vegetation index; WET, wetness; BSI, bare soil index; SI, salinity index; RSEI, remote sensing-based ecological 
index. 


Table 5 Socio-economic data of the Hotan Oasis 


Total Total power of Gross 


; . Cultivated Grain : Primary Secondary Tertiary 

Year Paonia ei land area production ded en industry GRP industry GRP industry GRP 

persons) (10,000 kW) (1000 hm?) (t) (10,000 CNY) (10,000 CNY) (10,000 CNY) (10,000 CNY) 
1992 90.10 17.50 81.84 387,358 81,029 50,761 11,648 18,620 
1995 96.60 18.40 82.64 404,569 160,780 105,131 21,612 34,035 
2000 110.00 23.00 91.17 509,114 178,411 101,566 22,614 54,231 
2005 120.00 26.00 91.41 579,342 321,140 134,080 75,448 111,612 
2010 135.00 33.00 97.11 674,272 696,522 210,118 141,565 344,839 
2015 155.00 54.00 - - 1,612,653 403,842 266,778 942,033 
2019 171.28 - 123.25 - 2,630,573 451,317 422,743 1,756,513 


Note: Kunyu City was established in 2016 and the data is incomplete. GRP, Gross Regional Product. The GRP is based on the prices of 


that year; - means that no statistics are available. 


areas have increased the size of impervious surfaces across the oasis. They have also contributed 
to the fragmentation of RSEI distribution. Results by Xue et al. (2021) indicated that 
anthropogenic factors as the main driving force led to the increased fragmentation in the 
landscape ecological pattern of the Hotan Oasis. Anthropogenic influences have also led to the 
expansion of the oasis, which is similar to the results of this study. 

From the macroscopic perspective of satellite remote sensing monitoring, we can assume that in 
the oasis ecosystem, the greater the greenness index, i.e., the greater the vegetation cover, the 
greater the surface humidity and the smaller the scale of bare ground and salinity, the higher the 
biodiversity and ecosystem complexity and the more excellent the eco-environmental quality. 
Moreover, the application of medium resolution remote sensing images is inevitably be affected by 
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mixed pixels. Therefore, more field validation data in some key areas in further work should be 
supplemented. 


5 Conclusions 


In this study, we used PCA to couple greenness, humidity, dryness and salinity on the GEE 
platform to construct a macroscopic and efficient monitoring method for eco-environmental 
quality based on improved RSEI. The spatial and temporal distribution of eco-environmental 
quality in the Hotan Oasis were quantitatively observed and its cause was analyzed in conjunction 
with LUCC. Results showed that the improved RSEI can reflect eco-environmental quality 
effectively of the oasis in arid area with relatively high applicability. NDVI and WET had a 
positive effect on eco-environmental quality of the Hotan Oasis. Meanwhile, BSI and SI had 
negative effects. The spatial distribution of eco-environmental quality for the Hotan Oasis is high 
in the middle, and low in the surrounding areas, decreasing from the middle to the surrounding 
areas of the oasis. Moreover, RSEI at the periphery of the oasis changes from low to high with 
increasing years, while the fragmentation of RSEI distribution within the oasis increases. This is 
mainly due to the impact of anthropogenic, including the expansion of artificial oasis into the 
desert, the replacement of desert ecosystems by farmland ecosystems, and the increase in the 
distribution of impervious surfaces. In arid areas with limited water and fertilizer resources, 
changing land use methods from quantitative expansion to high-quality development can reduce 
the risks that the oasis ecological environment may face. 
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